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Abstract 



Following the recent work of Caves, Fuchs, and Rungta [Found, of Phys. 
Lett. 14 (2001) 199], we discuss some entanglement properties of two-rebits 
systems. We pay particular attention to the relationship between entangle- 
ment and purity. In particular, we determine (i) the probability densities for 
finding pure and mixed states with a given amount of entanglement, and (ii) 
the mean entanglement of two-rebits states as a function of the participation 
ratio. 
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I. INTRODUCTION 



It has been recently pointed out by Caves, Fuchs, and Rungta H that real quantum 
mechanics (that is, quantum mechanics defined over real vector spaces [@-|5[]) provides an 
interesting foil theory whose study may shed some light on which aspects of quantum en- 
tanglement are unique to standard quantum theory, and which ones are more generic over 
other physical theories endowed with the phenomenon of entanglement. 

Nowadays there is general consensus on the fact that the phenomenon of entanglement 
is one of the most fundamental and non-classical features exhibited by quantum systems 
f|. Quantum entanglement is the basic resource required to implement several of the most 
important processes studied by quantum information theory fi-|H|, such as quantum telepor- 



tation [12| , and superdense coding |L3| . A state of a composite quantum system constituted 
by subsystems A and B is called "entangled" if it can not be represented as a convex lin- 
ear combination of product states. In other words, the density matrix p AB represents an 
entangled state if it can not be expressed as 

p ab = Y,Pkp£®pl (!) 

k 

with < pk < 1 and J2kPk = 1- On the contrary, states of the form ([I]) are called 
separable. The above definition is physically meaningful because entangled states (unlike 
separable states) cannot be prepared locally by acting on each subsystem individually |14| . 



The entanglement of formation provides a natural quantitative measure of entanglement 
with a clear physical motivation |T5|JT6[] . 

In standard quantum mechanics the simplest systems exhibiting the phenomenon of 
entanglement are two-qubits systems. They play a fundamental role in Quantum Information 
Theory. It should be stressed that the concomitant space of (mixed) two-qubits states is 
15-dimensional and its properties are not trivial. An explicit expression for the entanglement 



of formation of a two-qubits state p has been found by Wootters ||16|| . Wootters' celebrated 
formula has allowed for a systematic survey of the entanglement properties of the space of 
two-qubits states ||i7|-p0[ . 



Within quantum mechanics defined over real vector spaces, the most basic kind of com- 
posite systems are two-rebits systems. Rebits are systems whose (pure) states are described 
by normalized vectors in a two dimensional real vector space. A rebit may be regarded as 
the simplest possible quantum object ||. An explicit expression for the entanglement of 
formation of arbitrary states of two-rebits has been obtained by Caves, Fuchs and Rungta 

The aim of the present work is to explore numerically the entanglement properties of 
two-rebits systems. We pay particular attention to the relationship between entanglement 
and purity. 

The paper is organized as follows. In section II we review the CFR expression for the 
entanglement of formation of arbitrary two-rebits states, and discuss some of its immediate 
consequences. The relationships, for two-rebits systems, between the amount of entangle- 
ment and the degree of mixture are investigated in sections III. Finally, some conclusions 
are drawn in section IV. 

II. THE CFR FORMULA AND SOME OF ITS CONSEQUENCES 

Caves, Fuchs, and Rungta formula for the entanglement of formation of a two-rebits 
state p reads []J 

m = h (i±4E^) , (2) 

where 

h(x) = — x\og 2 x — (1 — x) log 2 (l — x), (3) 

and the concurrence C is given by 

C[p] = | tr(r) | = | tr(p a y ®<r v )\. (4) 

The above expression has to be evaluated by recourse to the matrix elements of p computed 
with respect to the product basis, =\ i) \ j), i,j = 0, 1. 
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We are also going to need a quantitative measure of mixedness. There are several mea- 
sures of the degree of mixture that can be useful within the present context. The von 
Neumann measure 

S = -Tr(p]np), (5) 

is important because of its relationship with the thermodynamic entropy. On the other 
hand, the so called participation ratio, 

is particularly convenient for calculations ||T7| , |l9f . 

A remarkable property of two-rebits states, which transpires immediately from the CFR 
expressions is that their square concurrence (and, consequently, their entanglement 

of formation) are completely determined by the expectation value of one single observable, 
namely, o y <S> a y . On the contrary, it has been recently proved that there is no observable 
(not even for pure states) whose sole expectation value constitutes enough information to 



determine the entanglement of a two-qubits state |22[. The operator a y ® o y has eigenvalues 
1 and —1, both two-fold degenerated. Let us denote by | 0i, 2 ) the pair of eigenvectors with 
eigenvalue 1, and | 03,4) the eigenvectors with eigenvalue —1, so that 

2 4 

cr y a y = ^2 \ 4>i) (4>i \ (pi) ((pi I • (7) 

i=l i=3 



A notable consequence of the CFR expressions (@-(|) is that there are mixed states of two 
rebits with maximum entanglement (that is, with C 2 = 1). For instance all states of the 
form 

P = p\ 0i)(0i I +{l-p)\<h)(<h I, (8) 

with < p < 1, are maximally entangled. Hence, for any participation rate within the 
range 1 < R < 2 there exist two-rebits states with maximum entanglement. We shall 
return to this point later, when we discuss the distribution of general two-rebits states in 
the (R, C 2 )-plane. 



III. ENTANGLEMENT VS. PURITY FOR ARBITRARY TWO-REBITS STATES. 



A. Measure on the Two-Rebits State Space 

In order to explore numerically the properties of arbitrary two-rebits states, it is necessary 
to introduce an appropriate measure p on the space Sr of general two-rebits states. Such a 
measure is needed to compute volumes within the space Sr, as well as to determine what is to 
be understood by a uniform distribution of states on Sr. In order to find a natural measure 
on S we are going to follow a line of reasoning akin to the one pursued by Zyczkowski et al. 
T7l , |T8|l in the case of two-qubits states. 



An arbitrary (pure and mixed) state p of a (real) quantum system described by an 
iV-dimensional real Hilbert space can always be expressed as the product of three matrices, 

p = RD[{\i}}R T . (9) 

Here R is an N x iV orthogonal matrix and -D[{Aj}] is an iV x iV diagonal matrix whose 
diagonal elements are {Ai, . . . , A^}, with 1 > Aj > 0, and J2i At = 1. The group of orthogonal 



matrices O(N) is endowed with a unique, uniform measure v [|2T|. On the other hand, the 
simplex A, consisting of all the real iV-uples {Ai, . . . , Aa?} appearing in is a subset of a 
(N — l)-dimensional hyperplane of 1Z . Consequently, the standard normalized Lebesgue 
measure Cn-i on IZ 1 ^^ 1 provides a natural measure for A. The aforementioned measures on 
O(N) and A lead then to a natural measure p on the set Sr of all the states of our (real) 
quantum system, namely, 

p = v x £jv-i- (10) 

We are going to consider the set of states of a two-rebits system. Consequently, our 
system will have iV = 4. All our present considerations are based on the assumption that 
the uniform distribution of states of a two-rebit system is the one determined by the measure 
flnp. Thus, in our numerical computations we are going to randomly generate states of a 
two-qubits system according to the measure C[3~0|) . 



B. Entanglement and Degree of Mixture. 



The relationship between the amount of entanglement and the purity of quantum states 
of composite systems has been recently discussed in the literature |P7-pO|. The amount 
of entanglement and the purity of quantum states of composite systems exhibit a dualistic 
relationship. As the degree of mixture increases, quantum states tend to have a smaller 
amount of entanglement. In the case of two-qubits systems, states with a large enough 
degree of mixture are always separable [T7J . To study the relationship between entanglement 
and mixture we need quantitative measures for these two quantities. As already mentioned, 
the entanglement of formation provides a natural quantitative measure of entanglement with 
a clear physical motivation [|l^,|l^] . 

The continuous line in Fig. 1 depicts the behavior of the mean entanglement of formation 
(E) of real density matrices, as given by the CFR formula, as a function of the participation 
ratio R. The dashed line in Fig. 1 corresponds to the mean entanglement of the same 
matrices, as given by Wootters' formula. In other words, in Fig. 1 the continuous line 
describes the mean entanglement of formation of the real density matrices when regarded as 
defined on a real vector space, while the dashed line describes the entanglement of formation 
of these same matrices when they are considered in the context of the standard complex 
vector space. We see that the CFR formula always gives, for the mean entanglement of 
formation, a value larger than the one obtained by recourse of the Wootters' expression. In 
this respect, our numerical results are fully consistent with the general arguments provided 
in!- 

The continuous line in Fig. 2 illustrates the behavior of the mean entanglement of 
formation (E) of real density matrices (given by the CFR expression) as a function of the 
participation ratio R. The dashed line in Fig. 2 shows the behavior of the mean entanglement 
of formation (E) of complex density matrices (given by Wootters' formula) as a function of 
the participation ratio R. 

The largest eigenvalue X m of the density matrix constitutes a legitimate measure of 
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mixture, in the sense that states with larger values of A m can be regarded as less mixed. 
Its extreme values correspond to (i) pure states (with A m = 1) and (ii) the density matrix 
jl (with X m = 1/4). In Fig. 3 we depict the mean entanglement (E) of all the two-rebits 
states with a given value of their maximum eigenvalue A m , as a function of this last quantity. 
The upper line corresponds to the CFR expression and the lower line to Wootters formula. 
Notice that, in the case of Wootters formula, the mean entanglement vanishes for A m < 1/3. 

We have also computed numerically the probability P(E) of finding a two-rebits state 
endowed with an amount of entanglement E. In Fig. 4 we compare (i) the distributions 
associated with two-rebits states with (ii) the distributions associated with two-qubits states 



which were recently obtained by Zyczkowski et al. ||T7| , |23|| . Fig. 4a depicts the probability 
P(E) of finding two-qubits states endowed with a given entanglement E (as computed with 
Wootters' expression). The solid line correspond to arbitrary states and the dashed line to 
pure states. In a similar way, Fig. 4b exhibits a plot of the probability P(E) of finding two- 
rebits states endowed with a given entanglement E (as computed with the CFR formula). 
The solid line correspond to arbitrary states and the dashed line to pure states. Comparing 
Figs 4a and 4b we find that the distributions P(E) describing arbitrary states (that is, both 
pure and mixed states) exhibit the same qualitative shape for both two-qubits and two- 
rebits states: in the two cases the distribution P(E) is a decreasing function of E. On the 
contrary, the distribution P(E) corresponding to pure two-rebits states differs considerably 
from the one associated with pure two-qubits states. The probability distribution P(E) for 
pure states of two-rebits reaches its maximum value for separable states (E — 0), and it is a 
monotonous decreasing function of the entanglement of formation E. On the contrary, the 
distribution corresponding to pure states of two-qubits is an increasing function of E for low 
values of the entanglement, and decreases with E for large enough values of this variable. It 
adopts its maximum value for an intermediate value of E. The general conclusion that we 
may draw from Fig. 4 is that the two curves representing the distributions P(E) associated 
with (i) pure states and (ii) arbitrary states do not differ, in the case of two-rebits states, as 
much as they do in the case of two-qubits states. 



The distribution P(E) for pure two-rebits states can be obtained analytically. Let us 
write a pure two-rebits state in the form 

4 



i=l 



(11) 



where 



E 

4=1 



i, a g ?e. 



(12) 



The states (| 0«), z = 1, . . . , 4) are the eigenstates of the operator a y ®a y , in the same order 
as in equation (|7|). The four real numbers q constitute the coordinates of a point lying in the 
three dimensional unitary hyper-sphere S3 (which is embedded in TZ 4 ). We now introduce 
on 53 three angular coordinates, 0i, 02, and 9, defined by 

C\ = COS 9 COS 0i, 

c 2 = cos 9 sin 0! , 
c 3 = sin ^ cos 02, 
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C4 = sin 6 1 sin 02, < 9 < — , < 0i,02 < 27r. 



(13) 



In terms of the above angular coordinates, the concurrence of the pure state | *f>) is given by 



C = I (a y ® cr y ) I = I cos 26* I 



(14) 



The element of volume on the three dimensional hyper-sphere is sin 9 cos 9d9d<p\d<p2- Thus, 
the total volume associated with a small interval d9 is 



dv = An sin 9 cos 9 d9. 



(15) 



Inspection of equations (0) and (|T5| ) allows one to deduce that the probability density P(C) 
of finding a pure two-rebits state with concurrence C is 



P(C) 



1 
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dv 



d9 



dO 



dC 



1. 



(16) 



The concomitant probability density P(E) of finding a pure state with entanglement of 
formation E is then equal to 



S 



(17) 

where dE/dC is to be computed from expressions (0) and (|J). It can be verified from 
equation ( |ITD that the limit value of P(E) associated with states of maximum entanglement 
is P(E = 1) = In 2. This value corresponds to the horizontal line in Figure 4b. 



P(E) 



dE 



dC 



C. Maximum Entanglement Compatible with a Given Degree of Mixture. 

We are now going to determine which is the maximum entanglement E m that a two- 
rebits state with a given participation radio R may have. Since E is a monotonic increasing 
function of the concurrence C, we shall find the maximum value of C compatible with a 
given value of R. In order to solve the concomitant variational problem (and bearing in 
mind that C = \ (cr y (g) a y ) | ), let us first find the state that extremizes Tr(p 2 ) under the 
constraints associated with a given value of (a y <g> a y ), and the normalization of p. This 
variational problem can be cast as 
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Tr(p 2 ) + /3(a y <g> a v ) - aTr(p) = 0, (18) 



where a and (3 are appropriate Lagrange multipliers. The solution of the above variational 
equation is given by the density matrix 

1 



Pr, 



al - (3(a y <g> (j y ) . (19) 



The value of the Lagrange multiplier a is immediately determined by the normalization 
requirement, 

Tr(p) = l =► a=i (20) 

Consequently, those two-rebits states yielding the extremum values of Tr(p 2 ) (and also the 
extremum values of R = 1/Tr(p 2 )) compatible both with normalization and a given value 
of (<j y <g> <j y ) are described by the density matrix 
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Pm = \l ~ 20(°tf® °v), ( 21 ) 



with the Lagrange multiplier (3 lying in the interval 



1 1 

2' 2 



(22) 



Density matrices of the form (pTj) , corresponding to negative values of /3, have been consid- 
ered in [p]], although not in connection with the variational problem that we are discussing 
here. 

In terms of the parameter (3, the expectation value (a y ® a y ), the concurrence squared 
C 2 , and the participation ratio of the statistical operator p m are given by 

(cr y ® a y ) = -2(3, 



C 2 = {a y ®a y ) 2 = 4(3 2 , 

R = t£p> ^ 

Hence, the maximum value of R compatible with a given value of C 2 is given by 



Rm{c 2 ) = y^c 2 ' ( 24 ) 



Rm{C 2 ) is a monotonic decreasing function of C 2 and adopts its values in the interval 
2 < R < 4. This implies that, within this range of i?-values, the maximum value of C 2 
compatible with a given value of R is the one obtained when solving Eq. (p4[) for C 2 , 
namely, 

C 2 = i - 1. (25) 

On the other hand, for 1 < R < 2 there always exist density matrices of maximum entan- 
glement (that is, with C 2 = 1). As a consequence, the maximum value of C 2 compatible 
with a given value of R is given by 

In Fig. 5 we plot (in the (R, C 2 )-plane) one million numerically generated random two- 
rebits states. The solid line corresponds to the maximum concurrence squared C^, for a 
given value of the participation radio R, as given by Eq. (|26|) . 
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IV. CONCLUSIONS 



We have explored numerically the entanglement properties of two-rebits systems. We 
paid particular attention to the relationship between entanglement and purity. We have 
computed numerically the mean entanglement of formation of two-rebits systems (as de- 
termined by the CFR formula (|2J)) as a function of the participation ratio R. We have 
also determined numerically the probability densities P(E) of finding (i) pure two-rebits 
states and (ii) arbitrary two-rebits states, endowed with a given amount of entanglement E. 
Furthermore, we surveyed the distribution of general two-rebits states in the (R, C 2 )-plane. 
In particular, we determined analytically the maximum possible value of the concurrence 
squared C 2 of two-rebits states compatible with a given value of the participation ratio. An 
interesting feature that deserves special mention is that, with regards to the probability of 
finding states with a given amount of entanglement, the difference between mixed states and 
pure states is much larger for qubits than for rebits. 



It would be interesting to perform, for quaternionic quantum mechanics |p4}| , an analysis 
similar as the one done here. 
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FIGURE CAPTIONS 



Fig. 1- Mean entanglement of formation (E) of real density matrices as a function of the 
participation ratio R. The continuous line corresponds to the CFR formula. The dashed line 
corresponds to the mean entanglement of the same matrices, as given by Wootters' formula. 

Fig. 2- The continuous line shows the behavior of the mean entanglement of formation (E) 
of real density matrices (given by the CFR expression) as a function of the participation 
ratio R. The dashed line shows the behavior of the mean entanglement of formation (E) 
of complex density matrices (given by Wootters formula) as a function of the participation 
ratio R. 

Fig. 3- Mean entanglement (E) of all the two-rebits states with a given value of their 
maximum eigenvalue X m , as a function of this last quantity. The upper line corresponds 
to the CFR expression and the lower line to Wootters formula. Notice that, in the case of 
Wootters formula, the mean entanglement vanishes for \ m < 1/3. 

Fig. 4a- Plot of the probability P(E) of finding two-qubits states endowed with a given 
entanglement E. The solid line correspond to arbitrary states and the dashed line to pure 
states. 

Fig. 4b- Plot of the probability P(E) of finding two-rebits states endowed with a given 
entanglement E. The solid line correspond to arbitrary states and the dashed line to pure 
states. The horizontal line corresponds to the limit value P(E = 1) = In 2 of the probability 
density associated with pure two-rebits states. 

Fig. 5- Plot of in the (R, i?)-plane of one million random numerically generated two-rebits 
states. The solid line corresponds to the maximum entanglement E m , for a given value of 
the participation radio R. 
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fig. 2 




fig. 4 




